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Abstract 

Evolution of genus density is calculated from Gaussian initial conditions using Zel'dovich ap- 
proximation. A new approach is introduced which formulates the desired quantity in a rotationally 
invariant manner. It is shown that normalized genus density does not depend on the initial spectral 
shape but is a function of the fluctuation amplitude only. 
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1 INTRODUCTION 



In the near future, we will obtain enormous information of the large-scale structure (LSS) of 
the Universe from several red-shift surveys. The ESO Slice Project (Vettolani et al. 1995) and 
the Las Campanas Red-Shift Survey (Shectman et al. 1995), have recently been completed, and 
data analysis of them are now being carried out energetically. The Sloan Digital Sky Survey 
(SDSS)(Loveday 1996) to be completed in the beginning of the next century will provide us with a 
three dimensional map of the Universe with about one million galaxies, which means that the data 
available will be one hundred times bigger than those present. 

Before the arrival of these data, it is important to furnish statistical measures to analyze them 
so that we will be able to probe the seed of density fluctuations as well as their evolution. A 
number of measures have been proposed to characterize statistical properties of LSS. The most 
commonly used quantity is the two-point correlation function (Totsuji & Kihara 1969), with which 
we can find the amplitude of fluctuations on various scales but not their statistical distributions. 
The counts-in-cells statistic, including the void probability (White 1979), directly measures the 
statistical distribution of galaxies and various theoretical models based on physical or mathematical 
arguments have been proposed to fit the observational data. Although this measure contains 
mathematically full information of the statistics in principle, it is difficult to relate the count 
analysis with the visual image or connectivity of galaxy clustering such as filamentary networks, 
sheet-like or bubble- like structures, etc.. As a statistical measure to characterize such a topological 
structure of galaxy distribution, the genus number has been used in the analysis of recent redshift 
surveys (Gott, Melott, & Dickinson 1986; Gott, Weinberg, & Melott 1987; Weinberg, Gott, & 
Melott 1987; Melott, Weinberg, & Gott 1988; Gott et al. 1989; Park & Gott 1991; Park, Gott, k 
da Costa 1992; Weinberg & Cole 1992; Moore et al. 1992; Vogeley, Park, Geller, Huchra, & Gott 
1994; Rhoads, Gott & Postman 1994). 

Theoretically, however, the value of the genus had been calculated only for the random Gaussian 
field (Adler 1981; Doroshkevich 1970; Bardeen Bond, Kaiser, &; Szalay 1986 (hereafter BBKS); 
Hamilton et al. 1986) for a long time except for some restricted cases (Hamilton 1988, Okun 1990). 
Although the random Gaussian distribution is commonly assumed as an initial condition, nonlinear 
gravitational evolution of fluctuations certainly changes the statistical distribution, and hence the 
Gaussian formula is not useful to trace evolution. The situation was somewhat improved recently 
as lowest-order correction to the Gaussian genus number was analytically obtained by one of us 
using the multi-dimensional Edgeworth expansion around the Gaussian distribution (Matsubara 
1994). Detailed comparison has also been done with the results of A^-body simulations, and it 
has been shown that the new formula fits the numerical data well in the semi-linear regime but 
not in the nonlinear regime (Matsubara &; Suto 1996). This is in accord with the fact that the 
one-point probability distribution function (PDF) based on the Edgeworth series no longer fits 
the counts-in-cells when the root-mean-square (RMS) value of the density contrast becomes as 
large as ~ 1/4 (Juszkiewicz et al. 1995; Ueda k, Yokoyama 1995). Thus it was desired to find 
an analytic expression of the genus number for realistic non-Gaussian distributions just as various 
non-Gaussian models have been proposed for counts-in-cells statistics. 

This was achieved by Matsubara & Yokoyama (1996, hereafter MY) who calculated the genus 
density in the case density field is characterized by a number of Gaussian variables such as the 
lognormal and the chi-square distributions. They have found that the lognormal formula fits 
evolution of a low-density cold-dark-matter (CDM) universe with a cosmological constant well. 
Although its agreement is excellent and the lognormal distribution is well motivated in the sense 
that it can be obtained from the continuity equation in the nonlinear regime but with linear or 
Gaussian velocity fluctuations (Coles &; Jones 1991), one cannot explain why it flts some models 
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but not others, since it is merely a fitting formula. 

In the present paper we calculate nonlinear evolution of genus statistics using the Zcl'dovich 
(1970) approximation (ZA) in order to clarify the dependence of the results on the initial spectral 
shape of Gaussian fluctuations. ZA describes time evolution of density field more accurately than 
the Eulerian perturbation theory in the semi-nonlinear regime. In particular, it has been claimed 
that ZA or its truncated version fit the pattern of the particles in iV-body simulations well (Coles 
et al. 1993). 

We make use of the Lagrangian description of ZA, in which time evolution of density field is 
easily described using initial fluctuations of gravitational potential, as well as the topologically 
invariant nature of the genus number which does not depend on coordinate choice. For further 
calculative simplicity, we introduce a new idea to deal with the time evolution of the genus number 
density. That is, we focus on the events of topology change of isovalue surfaces in 3-|-l-dimensional 
spacetime and relate it with the desired quantity with the help of the Poincare-Hopf theorem in 
differential topology. Thus we can formulate the genus density in a rotationally invariant manner 
and explicit evaluation becomes possible even if we must deal with 31 Gaussian variables in ZA as 
explained below. 

The plan of the present article is as follows. In §2 the basic concepts of genus statistics are 
reviewed and we give an expression when the density contrast is random Gaussian-distributed. In 
§3 a brief description of ZA is given. Then an important relation between the Lagrangian mapping 
and genus number is explained. In §4 we discuss a new formulation of time evolution of genus 
number density. We can easily understand that this formalism is a natural extension of a well 
known peak formalism in a three dimensional case. Numerical results are described in §5. §6 is 
devoted to conclusion and discussions. We add four appendices. In Appendix A, we present a 
mathematical foundation of our new formalism. More intuitive discussion is given in Appendix 
B. In Appendix C, calculative procedure is explained in detail. Finally in Appendix D, various 
perturbative formulae are given for comparison. 



2 GENUS CURVE IN THE GAUSSIAN DISTRIBU- 
TION 

Genus is a topological quantity defined by the number of the homotopy classes of closed curves 
that may be drawn on a surface without cutting it into two pieces. The surface we consider in 
this article is the isodensity contour of matter density field. We deal with the comoving number 
density of the genus on these surfaces. According to Morse theory there are a number of equivalent 
expressions of the genus density. 



r , . 1 [# (holes) — # (isolated regions)! 

genus density = '- ^ (1) 

volume 

— [^(maxima) -|- ^(minima) — #(saddle point)] 



2 volume 

For example, in the case of a 2-sphere, number of hole is 0, number of isolated region is 1, and 
genus number is —1 using formula (1). It has one maximum and one minimum along any direction, 
so genus number calculated from expression (2) is also —1. The genus number of a torus is 0, which 
we can easily confirm by the both formulae. Usually we use the second formula, and evaluate the 
number density of the stationary points to the direction of the third spatial axis, xs, together with 
the sign of the determinant of the Hesse matrix to assign the proper signature to each stationary 
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point depending on whether it is an extremum or a saddle. Thus the genus number density, G'(z^, t), 
of the contour S{x,t) = iycr{t), where 5 is density contrast and a is its RMS, can be expressed as 
follows (Doroshkevich 1970; Adler 1981; BBKS 1986). 

G{u,t) = -^(<5drc Wx,t) -iva(t))5drc(??l)<Jdrc(??2)|%|(CllC22 -C?2)), (3) 

where r/j(x, t) = dxi6{x, t), Oj(x, t) = dxidxjS{x, t), and (5drc(') is Dirac's delta function. Clearly this 
formula (^) lacks 0(3) symmetry due to the special treatment of the third axis. We will overcome 
this point later. In the language of differential topology, the Morse vector field is identified with 
the vector field {dx^x^, 8x2X3) in the above expression (|3[), where X3 is regarded as a function of xi 
and X2 implicitly defined by 5{xi, X2, X3{xi, X2);t) = i'a{t). (See Appendix A for the definition 
of mathematical terminology.) 

If 6 is isotropic random Gaussian, then we have simple relations between the variance of 5, rj, 
and C- (BBKS 1986). 

{6\x,t)) 
{r]i{x,t)r]j{x,t)) 

where crfit) and aK^) are defined, respectively, by 

al{t)^{[V6{x,t)f), and al{t) ^ {[V'6{x,t)f) . (5) 

Then we can easily calculate the statistical average in equation (|3|) to find (Doroshkevich 1970; 
Adler 1981; Gott, Melott, & Dickinson 1986) 

2 3 

We should notice that o"! has no influence on equation (P). This is because we only need sign 
information of the second derivative. 

MY studied the genus number density in the case matter density field departs from random 
Gaussian. They first evaluated genus curve in the lognormal density distribution, whose statistical 
property is characterized by a monotonic function, i.e. logarithm, of a random Gaussian variable. 
They also extended this kind of analysis to the case density field is characterized by a number of 
Gaussian fields through a function as in the chi-square distribution. Their method is also applicable 
to the present case of ZA in which density field is essentially characterized by six independent 
random Gaussian variables and their spatial derivatives as explained in the next section. 

3 ZEL'DOVICH APPROXIMATION 

3.1 Basic Properties of ZA 

Zel'dovich Approximation is known as one of the most reliable approximation to describe the weakly 
nonlinear evolution of density fluctuation (Shandarin & Zel'dovich 1989). It is a Lagrangian theory 
giving the shift vector of matter density field in terms of fluctuations in the initial gravitational 



a\t), (5(x,t)??,(x,t)) =0, (5(x,t)Cy(x,t)) 
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-6ij, {rii{x,t)Cjk{x,t)) = 0, 

{^ij^kl + ^ik^jl + ^il^jk), 



(4) 
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potential, $o> as a function of the initial position q. The relation between Lagrangian coordinate, 
q, and Eulerian coordinate, x, is given by 



X = q - Z)(t)Vg$o 



(7) 



where D{t) is a linear growth factor of small density fluctuations. The functional shape of D(t) is 
determined by cosmological parameters. In this article D is sometimes used as a time coordinate. 
ZA is confirmed to be a good approximation to the stage when displacement D(t)'Vq^Q is not so 
small, or until the density fluctuations enter into a nonlinear regime (Coles et al. 1993). 

By the conservation of mass, the local density contrast in ZA, 5za, is given by the Jacobian of 
the Lagrangian to Eulerian mapping. 



5zAi^,t) + l 



d{xi,X2,X3) 



d{qi,q2,q3) 
\{1 - DX^){1 - D\2){1 - D\;)\-^ 

1 - Dtia + — {(tra)^ - tra^} - D^deta 



(8) 



where we have defined a 



and Aj's are the three eigenvalues of the tensor a. 



Except for caustics, where density diverges and ZA breaks down, the Lagrangian mapping 
is a homeomorphism. With the Lagrangian coordinate, ZA describes local density profile in a 
very simple manner. In order to calculate the genus number density of an isodensity contour 
SzAi'^^t) = i'cr{t), it is convenient to introduce a monotonic function of 6zAi^,t) which may also 
depend on time, F[5^^(x, t),t], and consider a contour F[6za, t] = F[i^a{t),t] to evaluate the genus 
number of that surface which is clearly equivalent with the contour (5^^(x, t) = uait). Choosing an 
appropriate shape of F greatly simplifies subsequent manipulations. Specifically we choose F as 



F[5zA{-^,t),t\ 



1 



m 
Ml 



D 



l + <5zA(x,t), 
\{l - DXi){l - DX2){1 - D\M 



(9) 



Using @j one can regard F as a function of either (x, t) or (q, t), which is denoted by F = f{x,t) 
and F = /(q, i), respectively. Since we would like to obtain genus density in the Eulerian space 
we consider the contour /(x, = c below. Following MY (eq. [18]) one can easily show that genus 
density of this contour at D{t), Q{c,D), is related with G{u,t) by 



G{v,t) = g{F[va{t),tlD{t)) 



(10) 



3.2 Variables in ZA 

The fact that ZA contains more information, i.e. three eigenvalues, than the Eulerian perturbation 
theory makes it more reliable. On the other hand, calculation of genus number becomes much more 
difficult. In fact we must deal with the following variables. 

= d,^d,^^,, (11) 

fiijk = dq^dq.dq^^O, (12) 
lijkl = dq^dq.dq^dq^^O- (13) 
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We assume initial density fluctuations and so $o are random Gaussian-distributed. Their correla- 
tion functions are then given by 

{(y-ijOi-ki) = ^so^iSijSki + SikSji + SiiSjk), (14) 

= so'd%, (15) 

{oiijPklm) = 0, (16) 

{aijlklmn) = -Sl^4S«™„, (^^) 

{f3ijkl3inm) = si^df-li^^, (18) 

{(^ijklbnno) = 0, (19) 

{lij kl^mnop) '^2 d^f.^^^^, (20) 

in) 

where ^ is an n-dimensional completely symmetric tensor normalized as 

d^.r = n\. (21) 



In order to orthogonalize a and 7, we define new variables Uijki as 



s2 "{^^■^'='> 20^2''"^"'^ ^ 

As can be easily checked {oiijUkimn) vanishes and two body correlation of (jJijki is given by 

.4 fiQ„ 2 



K^ijkl^mnopl — (^ijklmnop "'mnop{ij"l^l} ]^Qg2 "■ijkl^mnop- l^^j 

Since a, /?, and a; are all Gaussian distributed, we can express the probability distribution 
function (PDF), Pi^a, in terms of sq, si, and S2 defined above. Denoting a, (3, and u 

correctively by yi, ■ • ■, 2/31, where six a's, ten /3's, and fifteen w's make up 31 variables in total, the 
PDF reads 

PLiyi,---,y3i)dyi---dysi = p-^-^-^dy^ . . -dyai, (24) 

Q = lJ2yii^'%yj^ (25) 

^fij = {yiVj)- (26) 

Note that this is a PDF in the Lagrangian space. We must calculate the average using the Eulerian 
PDF, PE{yi; • • • , 2/31)7 which would give a volume fraction in the Eulerian space in a ergodic system. 
Pe and Pl are related by 

J(x- q) 

PE{yi, ■ ■ yzijdyi ■ ■ ■ dy^i = \ .. PUyi, • • • , y3i)dyi ■ ■ ■ dysi, (27) 

where J(x; q) = d{xi, X2, X3) /d{qi, q2, qs) is the Jacobian and (J(x;q)) is its average with respect 
to the Lagrangian PDF Pl (Kofman 1994). Since only the combination of Dsq has a physical 
meaning, we may set so = 1 to fix normalization of D. However, we will put so explicitly when 
appropriate below. Bernardeau &; Kofman (1995) derived an analytic formula of one-point density 
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PDF using ZA. Their result is rather comphcated due to the divergence of density fields at the 
caustics points. Furthermore RMS amplitude of fluctuation diverges in ZA. Hence we will adopt 
the variance in linear theory, CTi{t) = \/120so-D(t), to normalize the amplitude of 6. 

Since / reduces to a linear combination of a in the limit D — > 0, /(x, 0) = tra, it is also 
random Gaussian then. Therefore at D = 0, when the Eulerian coordinate coincides with the 
Lagrangian counterpart, the genus density, Q{c,0) is easily obtained in the form given in § 2, 

When D is finite, /(x, D) is now a nonlinear combination of a and calculation would become 
extremely complicated and tedious if we used the usual expression (|3|) without rotational symmetry 
due to the fact that we tried to count maximums, minimums, and saddle points of an isodensity 
surface along the j;3-direction. In the next section we introduce a new approach which greatly 
simplifies evaluation. 



4 3+1 FORMALISM 

4.1 Basic Idea 

Here we present a new formulation based on the fact that time evolution of genus number density 
can be expressed in an 0(3) invariant manner. We know genus number density Q{c,0) on the 
three dimensional spacelike hypersurface Sq at t — > or D{t) — > as in §2. We can trace its 
time evolution making use of the fact that the points (events) where genus number changes ( 
P, Q, and R, in Figure 1) correspond to a minimum, a maximum, or a saddle point, namely, 
stationary points of an isovalue hypersurface of /(x, D) = c along the time- or Z)-direction in four 
dimensional spacetime (xi, X2, x^, D). The number density of these stationary points can be 
statistically expressed as in three dimensional case following similar argument to derive equation 
To evaluate the statistical average of the frequency of these events is essential in our new 
formalism. This idea is mathematically supported by a relationship derived from Poincare-Hopf 
theorem, an important theorem in differential topology (see Appendix A): 

2-md(X, V) = x(Si)-x(So). (29) 

Here V is a three dimensional manifold which has two two-dimensional boundaries, Sq and Si. In 
our case V is a three dimensional isodensity contour embedded into the four dimensional spacetime, 
and Sq and Si are the initial and the final two-dimensional spatial section of constant time or D, 
respectively. ind{X, V) is index of Morse vector field X in V, and x(^j) is Euler number in each 
spatial section Sj, which is related to genus number g as x = ~2g. This formula shows that 
genus number in the final surface Si of interest is determined by that of initial surface Sq and the 
information of stationary points in three dimensional region V. 
In the present case, the Morse vector field X is taken as, 

X = {d^,D„ d^,D,, d^,D,), (30) 

where Dc is a function of {xi, X2, x^) implicitly defined by 

/(Xi, X2, X3, Dc{xi, X2, X3)) = c. (31) 

As explained in Appendix B, we can classify stationary points into four groups by the signature of 
the Hesse Matrix TCij = d'^Dc/dxidxj. 
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type 


signature 


A5 


1 creation 


+ + + 


-1 


2 merging 


+ - - 


-1 


3 annihilation 




+1 


4 split 


+ + - 


+1 



1 <-> 3 and 2 <-> 4 are related to each other by time reversal. Here one should notice that the 
signatures and signs of the Hesse Matrix are related to changes of genus numbers, Ag, as, 

Jl(signature) = sgn(det(?t:)) = -Ag. (32) 

So we have only to evaluate four dimensional correspondence of equation (P) . 

After simple consideration, we find the net number density, Nst, of these stationary points 
including sign information in four dimensional space as 

Nst = — ^drc ^dvcidxif)Sdj:cidx2f)^dj:cidxsf)dDfdet{dxidxjf). (33) 

This is a natural extension of (^ to a three dimensional hypersurface. So the genus density Q on 
the three dimensional space Si on which D = Di is given by 

gU = c, D,) = G{f = c, 0) + j^' {Nst)^p^ dD, (34) 
where the average (• • ■)af3uj is to be calculated with the Eulerian PDF (eq. [p7[|). 

4.2 Method for Calculation of (Nst) 

In order to calculate {Nst) we first rewrite equation ( |33| ) in terms of f{q,D) = /(q — DVq^o, D). 
From 

dl 

dxi 

dxidxj 
dD 

where f^i = df /dqi etc., we find 

Nst = -Sdrc (/(q, D) - c) <^drc(/,i)5drc(/,2)5drc(/,3)/,Ddet(/,,,)| J(x; (36) 
Hence we obtain 

{Nst) = J NstPE{y)d^^y 

- I Sdrc (/(q, D) - C) ^drc(/,l)^drc(/,2)(^drc(/,3)/,Ddet(/,,,)PL(j/)ci=^^?/ 

(■/(x;q)) ■ ^''^ 

The above equality can be interpreted as follows. The number of stationary points in a certain 
region remains invariant whichever coordinate, Lagrangian or Eulerian, one uses. So its density 
changes only by the volume factor represented by the Jacobian (J(x; q)), which has been calculated 
by Kofman et al. (1994). 



f.lm 



dqi 
dqi dqr, 



dxi dxj 



+ f,i 



d'^qi 
dxjdxi 



(35) 



df ^ d^odqm ^^ 

dD dqi dxi 
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In the above formulation 32-dimensional integration must be accomplished in principle including 
/^-integral. Making full use of the rotational symmetry, however, we can perform most of these 
integrals such as those over (5 and uj algebraically following the scheme introduced by MY. 

First we note that thanks to the 0(3) invariance Ngt does not depend on the Euler angles of 
the transformation to principal axes of aij, and we may assume that aij is diagonalized from the 
beginning without loss of generality. So among six integrals over Uij, those corresponding to the 
angular integrals are trivial (see Appendix B of BBKS). Consequently we can replace 

Oil! = Al, a22 = A2, 0:33 = A3, ai2 = Q!23 = "31 = 0. (38) 

We now turn to averaging over (3 and oj whose details are given in Appendix C. The de- 
terminant of f^ij in Nst has 3! = 6 terms. Thanks to the 0(3) symmetry, we may identify 

/,ii/.23/,32; /,22/,i3/,3i, and /,33/,i2/,2i with each other in the expression to be averaged. We 
can therefore rearrange det/^jj as follows. 

det/,jj- ^ /,ll/,22/,33 - 3/,ll/,23/,23 + 2/,l2/,23/,31- (39) 

Furthermore, since the Gaussian average of odd multiple is equal to zero, we can drop odd terms 
of ujijki or Pijt;. Using these facts it is not so hard to calculate the integrals. 

First, integrate over ujijki- We only have to replace ujuj by {ujuj), terms containing ujujuj or uj by 
0, and do not change the terms with no co. Next, integrate over (5ijk- Using three delta functions 
<^drc(/i)) we can independently eliminate integration over Pni, P222, and l3s3s- As the rest of Pijk- 
integral (seven dimensional) is ordinary Gauss integral, we can algebraically manage them by such 
a formula 

iPaPbPcPd) = M,,M,^ + M,,M,a + M^^Mbc, (40) 

where M^^ is the correlation matrix {(3aPb)- This matrix is 7 x 7 but constituted from 2 x 2 ® 2 X 
2©2x2©lxl blocks due to our choice of diagonalized a. 

Finally we numerically integrate over Aj and D. The PDF of Ai, A2, and A3, has been given by 
Doroshkevich (1970). 

P(Ai,A2,A3)dAidA2dA3 = J^|(Ai-A2)(A2-A3)(A3-Ai)| 



487rcro 
X exp 



cro^ V 2 



d\id\2d\z, (41) 



where we have defined 



Ji = A1 + A2 + A3, (42) 
J2 = A1A2 + A2A3 + A3A1, (43) 
(To = v^so- (44) 

There is a subtle point in the integration over Aj, which is the emergence of caustics where 1 — DXi 
vanishes. Since these points make both ZA and our homeomorphic picture break down, our result 
would not be trustworthy if the region 1 — DXi < should affect the final result. For this purpose 
we integrate over Aj with two different ways, one over the full range of Aj and the other over the 
region 1 — DXi > only. These two results should agree with each other for our analysis to be 
justified. Note that even though our homeomorphic picture breaks down at caustics we do not 
encounter any divergent quantities in our analysis there unlike the case of one-point PDF studied 
by Bernardeau & Kofman (1995). 
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Because of 5drc(/ — c) in equation (p8|), we can eliminate Ai. So the number of dimensions to 
integrate is only three, A2, A3, and D. The function to be integrated is complicated very much, 
whose abbreviated form is given at the end of Appendix C, but it is easily dealt with numerically 
for the rapid convergence due to the overall Gaussian-like factor from the PDF of A. 



5 RESULTS 

Before presenting the numerical results we comment on the functional dependence of the final 
results of 

g{f = c, D,) = Gif = c, 0) + J^""" {N,t)^p^ dD. (45) 

As pointed out before, equation (p5|) does not contain S2- This is because we only need sign 
information of second derivative and not their amplitude. Furthermore in ZA case, si appears in 
equation (^) only as an overall factor si^ / sq^ as in the Gaussian formula (eq. Q). This implies 
an interesting fact that for any value of si, every genus curve has a similar shape parameterized 
only by D since we have set sq = 1- Therefore the ratio ^(c, Di)/Q{Q, Di) depends on D and c 



only and not even on si. From equation (10) genus density of the surface 6 = vai{ti) is given by 



We now proceed to the numerical results carrying out numerical integration with different final 
time (i.e., 5^Di)- We can arbitrarily chose the normalization of so as pointed out before. On 
the other hand si changes only the scale of vertical coordinate for the overall factor. Since we 
concentrate on the ratio G{v., Di)/G{0, Di) we can also put si = 1 in the calculative process. 

In Figures 2(a)-(d), we plot the genus curves at ai = 0.1, 0.22, 0.5, and 1.0. For comparison, 
we have also plotted the results of two perturbative formulae, second-order Eulerian perturbation 
(eq. |0l|) and ZA to the same order (eq. |102| ) . These curves show that all of them are very close 
to each other as expected in the region where density contrast is small (c'"i|i^|< 1/4), but deviate 
considerably in both ends (i7/|j^|> 1/4). In the same figures we have also depicted contribution of 
the first term in equation (|45| ) to show the importance of time evolution of stationary points. 

In Figure 2(d) with fx^ = 1 we have depicted two curves of ZA result. The solid line corresponds 
to the case integration range of Aj has been restricted to the region 1 — DXi > 0, while short- 
dash-dotted line represents the case integration has been done over the full range of Aj. Their 
discrepancy is at most about 10% indicating the accuracy of our results. For the cases with smaller 
ai these two curves practically coincide with each other so only one of them has been depicted in 
Figures 2(a)-(c). 

Our results show clearly that the peak of the genus curve based on ZA shifts to right or to the 
region u > 0. The same tendency is observed in the result of Eulerian perturbation theory without 
smoothing (Matsubara 1994). However, the peak of the genus curve shifts to left in the latter 
theory once uniform smoothing such as a Gaussian smoothing with a specific width, say 8/i~^Mpc, 
is applied. The results of CDM type A^-body simulations also exhibits shift to left after the same 
smoothing. Unfortunately we cannot apply the same sort of smoothing to the final configuration in 
our result because we are working in the Lagrangian coordinates which coincides with the Eulerian 
counterpart only at the initial epoch and because smoothing and time evolution do not commute 
with each other. 

The information of initial condition must be contained more directly in Lagrangian description 
than in Eulerian description. So some of it might have been lost in Euler-type smoothing. The 
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discrepancy of the shift of the peak might support this view imphcitly. Therefore a new, more 
adequate smoothing method is desired. This issue is under investigation now and we do not go 
into it further here. 

6 CONCLUSION 

We studied the gravitational evolution of genus number density of isodensity surface in terms of 
ZA. While ZA is more reliable than the Eulerian perturbation theory to describe the density profile 
in semi-nonlinear regime, it is much more difficult to calculate genus density in the former which 
depends on more variables. We have overcome this difficulty by introducing a new approach in 
which we focus on evolution of genus number in the 3+1 dimensional spacetime rather than dealing 
with that at a specific time from the beginning. We can formulate the former in a rotationally- 
invariant manner so that its calculation is easier and relate it with the desired quantity making 
use of the Poincare-Hopf theorem in differential topology. Intuitively, this formula is based on the 
fact that events of topology change of equal-time spatial isodensity surfaces can be identified with 
stationary points of the corresponding hypersurfaces in four dimensional spacetime as explained in 
Appendix B. Without such a new idea we would have to have dealt with 31 dimensional numerical 
integration instead of just three. 

We also mention that the Lagrangian nature of ZA is very suitable for our formulation. On 
the other hand, that made it impossible to apply uniform smoothing to the final result and so 
we have been unable to compare our result with observations or simulations. In the Lagrangian 
description information on the initial condition is reflected more directly hence it is desirable to 
invent a new smoothing method of discrete data which is compatible with Lagrangian theories in 
order to find out the nature of primordial fiuctuations. We might obtain such a means by changing 
the smoothing width as a function of local density which is often adopted in smoothed particle 
simulations. However, we should clarify theoretical basis before drawing final conclusion. This 
issue is under investigation now. 

Apart from problems of smoothing we have shown that normalized genus density is a function 
of D only and does not depend on the initial spectral shape as long as initial fluctuation is Gaussian 
and ZA is applicable. 

NS would like to thank Professor H. Sato for his continuous encouragements. This work was 
partially supported by the Japanese Grant in Aid for Science Research Fund of the Ministry of 
Education, Science, Sports and Culture No. 08740202(JY) and No. 5391(MS). 

APPENDICES 
A POINCARE-HOPF THEOREM 

In this appendix we review a mathematical basis of our new formulation. The following Poincare- 
Hopf theorem (Milnor 1965) is essential for our analysis. 

Theorem I (Poincare-Hopf) Let M"' be a compact n-dimensional (n > 2) C°° manifold. X is any 
C°° vector field with at most a finite number of zeros, which are the points with vanishing 
magnitude of X, satisfying the following two conditions. 

(a) The zeros of X are contained in IntM^. 
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(b) X has outward directions on dAi^. 

Then the sum of the indices of X at all its zeros is equal to the Euler number x of A^", 

ind{X) = xiM"")- (47) 

Here the index of the vector field X at a zero, P, is defined as follows. Let Xa{x) be the components 
of X with respect to local coordinates {x"} in a neighborhood of P. Set Va{x) = Xa{x)/\X\. If we 
evaluate v on a small sphere centered at x{P), we can regard Va{S^^^) as a C°° mapping from S^~^ 
into S^~^. The mapping degree of this map is called the index of X at the zero P. For example, 
if the map is homeomorphic, the mapping degree of the orientation-preserving (reversing) map is 

+1 (-!)• 

In discussing genus statistics we identify the vector field X with a gradient of a C°° function h 
called the Morse function. If P is a non-degenerate critical point of this function with Morse index 
A, which is the number of the negative eigenvalues of Hesse matrix dadbh, then the index of the 
gradient field dah at P is (—1)'^. In the present paper Dc{xi,X2,xs) defined in equation ( [31] ) serves 
as the Morse function. 

We treat three dimensional manifold embedded in a four dimensional spacetime manifold as 
an isovalue contour. The three dimensional manifold has two two-dimensional boundaries as an 
initial boundary and a final boundary. For such a manifold, we use the following modification of 
the Poincare-Hopf theorem. Now we consider odd-dimensional manifold with two boundaries So,i- 

Theorem II (Sorkin 1986) Let M"" be a compact n-dimensional (n > 2 is an odd number) C°° 
manifold with Sq U Si = dAi"' and Sq H Ei = 0. X is any C°° vector field with at most a 
finite number of zeros, satisfying the following two conditions: 

(a) The zeros of X are contained in IntM"^. 

(b) X has inward directions at Sq and outward directions at Si. 

Then the sum of the indices of X at all its zeros is related to the Euler number of Sq and Si 
as 

x(Si) - x(So) = 2 • ind{X). (48) 

Proof: Let be a manifold obtained by gluing 7\A = Sq x [0, 1] and M.'^ at Sq x [1]. Also let X 
be an extension of X to which has at most a finite number of zeros and is outward directed on 
So X [0] C J\f . Since X and its restriction Xj^ to N are outward directed on the entire boundaries 
of and N , respectively, we can apply the theorem I to them. Then we have 

ind{X) = xiM^, (49) 

and 

ind{X^) = x{^f). (50) 

Since A^" and are diffeomorphic, xi-M.'^) = x(-^")- Using xi-^ x B) = xi-^)xiS) and 
ind{X) = ind{X) + ind{Xj\f), we obtain 

indiX) = xiMn-x{^o). (51) 

Similarly we find 

ind{-X) = x{Mn-x{^i), (52) 
by interchanging the roles of Sq and Si. Since ind{—X) = {—l)^ind{X) (inversion of 5*""^ reverses 
orientation only if n is odd), we have equality (^8|). 

Finally we consider time-vector field X = VDc- As mentioned above, the index of the vector 
field X is given by the number of negative eigenvalues of a corresponding Hesse matrix dadbDc- 
Therefore ind{X) in theorem II is evaluated in the same manner as in Doroshkevich (1970). 
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B STATIONARY POINTS IN 3+1 SPACETIME 



In this appendix we present an intuitive explanation of the relation between time evolution of 
genus number of two dimensional surfaces in three dimensional space and stationary points along 
the time axis on three dimensional hypersurfaces consisting of their trajectories in four dimensional 
spacetime. These points correspond to events where genus number of a two dimensional contour 
changes. They are classified into four groups by the signature of the Hesse matrix. The change of 
the genus number turns out to be minus of the sign of the determinant of that matrix as explained 
below. 

First we consider the case when a new isovalue surface is created at a spacetime point P with 
time t = tp, namely, when genus number changes by —1 as shown in equation (1). We introduce a 
local coordinate with its origin at P and principal axes denoted by X, Y, and Z in Figure 3-1. In the 
upper figures emergence of a new contour of f{X, Y,Z;D) = cor D = Dc{X, Y, Z) at tp is depicted, 
showing time evolution from tp — e to tp + e (e>0). Lower- left figures represent how intersecting 
points of the surface to each principal axis move with time, which is summarized in the lower-right 
figure showing evolution of the contour along each axis. In the present coordinate system we have 
only to investigate the sign of second derivatives d'^Dc/dX.'^, d'^Dc/dY'^, and d'^Dc/dT? of the three 
dimensional contour D = Dc{X, Y, Z) in four dimensional space. Apparently all of these derivative 
is positive, so the signature of the Hesse matrix is + + +, and its determinant is positive, too. 
Conversely when signature is + + + the type of the event is creation. Since annihilation is its time 
reversal the signature is . 

Next we consider the case when two surfaces merge at a spacetime point Q which is depicted 
in Figure 3-2 in which the origin is now taken at Q. Depending on whether the two surfaces are 
connected or disconnected this event is interpreted as merging of two isolated surfaces or creation 
of a hole. /,From the lower-right figure we can easily convince ourselves that the signature of Hesse 
matrix there is + H — , and genus number changes by +1. Since split, which corresponds to either 
increase of an isolated region or annihilation of a hole, is time reversal of this, the signature is 
H and genus changes by +1. 

We can summarize the above results that the change of genus number at each event is equal to 
minus the signature of the Hesse matrix. 



C CALCULATION SCHEME 



In this appendix, we briefiy explain the calculative procedure of {Nst)f^^ (see eq. |33|). This is 
essentially an extension of the scheme developed by MY. We perform uj- and /3-integral algebraically. 
The remaining integration over A and D has to be done numerically. Throughout this appendix 
we put So = 1. 

In the quantity to be averaged, 

Nst = -Sdrc (/(q, D) - C) 5drc(/,l)5drc(/,2)<5drc(/,3)/,Ddet(/,i,), (53) 

UJ appears only in the second derivative f^ij and (3 appears both in f^i and f^ij as 
df 

Pijk, (54) 
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J.ij o PiklPjmn + n 'Jijkl 

dakidomn oau 

— /-i o PiklPjmn ~l" ^ijkl 
daklOamn OUkl 



df (3 



jSi^a{ij^fc«} - -^si^ (an + 022 + 033) ((^ij^fci + ^ik^ji + ^i/f^jfe) j- • (55) 



daki [2 

Since aij has only diagonal components, Ojj = Xi5ij, we find (^jj = df /doij in a very simple form 
below. 

.n ^ ^ = ^ = -a-^A2)(l-m3)..„ (56) 

522 = -(l-i?Ai)(l-m3) =52, (57) 

533 = -{l-DXi){l-DX2)=g3, (58) 
9^J = (^/j). (59) 

Similarly hijki ^ d"^ f / duijOau reads, 

/iii22 = /12211 = -^1221 = -/i2112 = -0(1 - DAs) = /l3, (60) 
/11133 = ^3311 = -^1331 = -/13II3 = D{1 - DX2) = /l2, (61) 
^2233 = ^3322 = —^2332 = —^3223 = D{1 — DXi) = hi, (62) 

hijki = (others). 
As explained in the text, using 0(3) symmetry the determinant of f^ij can be replaced by 

det/,i,- => /,ll/,22/,33 - Knf% + 2/,12/,32/,13. (63) 

For notational simplicity we rewrite f3ijk as follows. 

fli = (3221, a2 = 13331, a3 = (3u2, a4 = /3332, (64) 
as = Pus, aQ = (3223, «7 = A23, (65) 

A = Pui, B = P222, C = P333- (66) 

We further introduce new variables and pi which are related with the first terms proportional 
to /3/3 and the third terms proportional to a in equation ([55|), respectively. They are defined as 
follows. 

ri = 2h3{Aai-al) + 2h2iAa2-al) + 2hi{aia2-a^), (67) 
r2 = 2/13(^03 - a?) + 2/i2(Sa4-a6) + 2/12(0304 -a?), (68) 

r3 = 2/ii (Cae - a|) + 2/i2(Ca5 - 02) + 2/13(0506 - 07), (69) 
3 

Pi = --s?(24Ai5i+4A252 + 4A353+4Ai52 + 4Ai53) 

+^s?(Ai + A2 + A3)(35i+52 + 53), (70) 
3 

P2 = - 2^i(24A252 + 4Ai5i + 4A353 + 4A251 + 4A253) 

+-s?(Ai + A2 + A3)(5i + 352 + 53), (71) 
5 
3 

P3 = --sf (24A353 + 4A252 + 4Ai5i + 4A35f2 + 4A351) 

+|s?(Ai + A2 + A3)(5i +92 + 353). (72) 
5 
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Then neglecting odd terms of to, which do not contribute to the average and using 0(3) symmetry 
further, we can replace /,ii/,22/,33 by 



/,ii/,22/,33 =^ rir2r3 + 3(ri + pi)gijUJij22gklWkm + 3pir2r3 + 3rip2P3 + PiP2P3- 
In the same manner, /,ii/23 and /,i2/,32/,i3 are written in simpler forms. 



/,ll/,23 =^ 
/,12/,32/,13 ^ 

where ri2, r23 and ri3 are defined as 



(ri + pi)gijUJij2-igki^km +P3'^i2?^i2, 
n2?'23n3, 



ri2 = /ii (01X4 - a2i? - 20507) + /i2(^a4 + 0203 - 20507) + /i3(Ai? - oioa), 

?'23 = /ll (-BC - O4O6) + /i2(a4a5 + - 2O2O7) + /l3(0306 + i?05 - 2O7O1), 

ri3 = /li(Coi + O6O2 - 2O7O4) + /l2(^C - O2O5) + /l3(A06 + O1O5 - 203O7). 



(73) 



(74) 
(75) 



(76) 
(77) 
(78) 



We perform uj- and /^-integral in order. In equation ( |73[) only the second term contains uj which 
is to be replaced by the corresponding variance (eq. ||2^). After some tedious algebra we obtain 
the following result. 



{gij^ij22gki^kmz) 



144s 



Hmgl + 1845152 + 1845193 + 92^253 + 138gi + 138^1). 



(79) 



Here we have dropped the terms proportional to S2, since this cancels with the similar term in ([7^) 
(see MY). 

We then move onto /^-integral which is of the form 

^ : exp --(3M-^f3] <5drc(/,i)5drc(/,2)<5drc(/,3) X (polynomial of (5). (80) 



(27r)VdetA^ 



Here ^ is the inverse of the correlation matrix of j3 whose components are given as follows. 



336s?7W"^ 
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04 
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05 


06 


07 
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-1/2 


3 


-1/2 
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-1/2 
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-1/2 


-1/2 
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-1/2 


3 
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-1/2 


-1/2 
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-1/2 


-1/2 
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-1/2 
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-1/2 
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-1/2 


-1/2 
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07 
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The factor 336 arises due to our normalization of the variance of 13. First we use three delta 
functions 5drc(/,i) and replace A, B, and C by 



^ = -{g2ai+ g3a.2)/gi, 

B = -(5103 +93a4)/5'2, 
C = -(5105 +52a6)/53, 



(81) 
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respectively. Then the remaining /3-integral is seven dimensional and we find 



(80) 



1 



VdetM 



91929?. (27r) 2 VdetTW J (27r) 2 VdetM 



: exp 



X (polynomial of a^). (82) 



Here Mj,- is a correlation matrix of Oj obtained from M. For example, for i,j = 1, 2 we find 



(aifli) 
(0102) 



192sf(9(72 + 35153 + 25! 



155i + 65152 + 65153 + 352 + 25253 + 353 



(«2ai) 



^ 96sf (35^ - 5152 - 35153 - 45253) 

155i + 65152 + 65153 + 35I + 25253 + 35I 
192sf(952 + 35152 + 25!) 



155i + 65152 + 65153 + 35I + 25253 + 35I 



(83) 
(84) 
(85) 



The other components are obtained similarly except for M77 = 48. As pointed out before this 
matrix is constituted from 2x2©2x2©2x2©lxl blocks. 

Finally the remaining Gaussian integral over is done by using the following relations of 
averages. 



{aiUjakai) 
{aiajakaittman) 



MijMklMmn + ■■■+ MinMjkMir. 



[15 terms). 



(86) 
(87) 



Thus we finish ui and /3 integration in equation ([73|). Other terms in equation ( [6^ ) can be dealt 
with in the same manner. Then we obtain the final form 



sti 



I3uj 



-(A1A2 + A2A3 + A3A1 — 2DX1X2X3) 



X V- 
X exp 



15 1 
3 (Ai + A2 + A3)^ + y (A1A2 + A2A3 + A3A1) I 



X 
X 



TT~^ ■ 2-75 . 3-15 . 5-1 . 7-3 
1 

120 

|(Ai-A2)(A2-A3)(A3-Ai)| 
(5i5253)"^<^drc(/ - c) VdetM 
|(nr2r3)^ + 3((ri)^+pi) 

+3(n)/3P2P3 +P1P2P3 
-3((ri)^+pi) sf (-2885? 

+2 (ri2r23r3i)^ 



144 



+ + ... 



(465? + ... + 13851) j +3pi (^2^3)^ 
25925I) - 3p3 (n2ri2)^ 



(88) 
(89) 



(90) 
(91) 



(92) 
(93) 

(94) 



Here the factor (^) comes from f^D term, and the factor (|96| ) from the PDF of A (Doroshkevich 
1970). the terms ( p^ ) through ( p^ ) arise from det/^jj. In the above {Nst)p^, we have already 

included the effect of (Euler) angular integral. So J^^^ J J J dDd?\{Nst) corresponds to the 
evolution of the genus number density. 



D PERTURBATIVE EVALUATION OF THE GENUS 
CURVE 

In this appendix, we present perturbative evaluation of the genus curve for comparison with our 
result. First we recall that Matsubara (1994) obtained correction to the Gaussian formula (eq. 
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Q) using the multi-dimensional Edgeworth expansion in the second-order Eulerian perturbation 
theory. The final result is 



1 



e 2 



(95) 



In the above expression, Hniv) = {—)^e'^^^'^{d/diy)^e ^'^ 1'^ is the n-th order Hermite polynomial, 
and S", T, and V are defined as 



5 = A(^')' 



(96) 



-^(V5-V5V25), 



respectively, which we call generalized skewness. They come from the three-point correlation func- 
tion that can be evaluated also by second-order perturbation theory (Matsubara 1994). The solution 
of the Eulerian second order perturbation theory of gravitational instability in Fourier space is, 



<^(2)(k,t) = Z)2| |^^(2vr)35L(pi +P2 -k)P2(Pl,P2,t)e(Pl)6(p2) 



where 



^2(Pi,P2,i) 



(1 + i^) + 



Pi • P2 V\ 



P1P2 \P2 Pi 



+ ^) +(1-K) 



Pi • P2 

P1P2 



(97) 



(98) 



In the above expressions, K(t) is a weak function of and A (Bouchet et al. 1992; Bernardeau et 
al. 1995; Matsubara 1995). In Einstein-de Sitter universe, K = 3/7. 
The parameters S, T and U are then evaluated as 



S = 2{2 + K), 



';{5 + 2K), U=-{3 + 2K). 



(99) 



The derivation of the above results of T and U is similar to the derivation of skewness, S (Bouchet 
et al. 1992; Peebles 1980). 

On the other hand, Grinstein &; Wise (1987) give an Eulerian representation of the Zel'dovich 
approximation. According to their result, the kernel P2 for Zel'dovich approximation is 



^'2(Pi,P2,i) 



1 k • pi k • p2 



2 p2 



P2 



1 + 



Pi • P2 /P2 , Pl\ , (Pl • P2 



— + — + 



^1^2 \Pl P2 



pIp2 



(100) 



Thus, K = Q corresponds to the Zel'dovich approximation in second-order Eulerian perturbation 
theory. 

The set of skewness parameters {S,T,U) is (34/7,82/21,54/35) for exact second-order per- 
turbation theory (of Einstein-de Sitter universe) and (4,10/3,6/5) for Zel'dovich approximation. 
Therefore, equation (p^) reduces to 
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1 — 1/ + av \ — I V V 
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21 



(101) 
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for exact second-order perturbation theory of Einstein-de Sitter universe, and 



2 \ 2 

e 2 



,7 5 n 

'53 



(102) 



for Zel'dovich approximation. 
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FIGURE CAPTIONS 



Figure 1 Schematic explanation of the evolution of the genus number in which one spatial di- 
mension is omitted. The genus number on the final surface Si is determined by the genus 
number on the initial surface Sq and information of stationary points in the spacetime in 
between. At point P a new region is created (in this picture the signature is ++). At point Q 

a region is annihilated ( ). At point R a region is split into two regions (— +)• In this three 

dimensional schematic picture, merger and split are not distinguished by signature (both have 

— h), but in the actual four dimensional case they are different; merger has H , while split 

has + H — . 

Figure 2 Normalized genus number density at ai = 0.1, 0.22, 0.50, and 1.0. Solid line corresponds 



to our new result of full ZA (eq. [45|), short-dashed line to second-order Eulerian perturbation 
formula G2ndi^) (eq. |101[ ), long-dashed line to second-order perturbation theory for ZA 
Gzei,2nd(f^) (eq. |102(| ), and short-dash-dotted line to the contribution of the first term of 
equation (|4^), namely, without the effect of time evolution of stationary points. In the last 
figure with ct/ = 1.0 long-dash-dotted line represents the full ZA result in which Aj-integration 
has been done up to infinity. 

Figure 3 Intuitive explanation for classification of stationary points in 3+1 spacetime. 
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Fig 2(a) 




Fig 2(b) 
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Fig 2(c) 




Fig 2(d) 
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